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Abstract. The proton's elastic electromagnetic form factors are calculated 
using an Ansatz for the nucleon's Poincare covariant Faddeev amplitude that 
only retains scalar diquark correlations. A spectator approximation is em- 
ployed for the current. On the domain of q 2 accessible in modern precision 
experiments these form factors are a sensitive probe of nonperturbative strong 
interaction dynamics. The ratio of Pauli and Dirac form factors can provide 
realistic constraints on models of the nucleon and thereby assist in developing 
an understanding of nucleon structure. 
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1 Introduction 

The pion's elastic electromagnetic form factor is accessible via a properly con- 
structed six-point quark Schwinger function, a fact exploited with modest success 
in lattice-QCD simulations, e.g., Refs. ^[^UHl- This Schwinger function is also the 
basis for continuum studies, among which those employing the Dyson-Schwinger 
equations (DSEs) [1103113 are efficacious [7j. The proton's form factors are acces- 
sible through an analogous eight-point quark Schwinger function, which is the 
starting point for lattice simulations, e.g. Refs. [HUH]. The fruitful extension of 
DSE methods to the calculation of this Schwinger function is a contemporary 
goal. As we will explain, a simple truncation that corresponds to a spectator 
approximation is currently in widespread use |10[ 11 2 j . Manifest covariance is 
a strength of this approach for it has long been apparent that in order to obtain 
an internally consistent understanding of proton form factor data at spacelike 
momentum transfers q 2 > M 2 , where M is the proton's mass, a Poincare covari- 
ant description of the scattering process is necessary |131 114] . This has recently 
been reemphasised in the context of constituent-quark models IT7|. 

The same interaction which describes the structure and properties of colour- 
singlet mesons also generates a quark-quark (diquark) correlation in the colour 
antitriplet (3) channel |18l I19j . Such correlations have recently been observed 
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in simulations of lattice-QCD 20 a . While diquarks do not survive as asymptotic 
states |2H 1221 155] : viz., they do not appear in the strong interaction spectrum, the 
existence of strong quark-quark correlations provides a foundation for viewing the 
nucleon as a quark-diquark composite. This picture can be realised via a Poincare 
covariant Faddeev equation [23], in which two quarks are always correlated as a 
3-diquark, and binding in the nucleon is effected by the iterated exchange of roles 
between the dormant and diquark-participant quarks, and through the action of 
a pion "cloud" j2S]. 

Upon solving the Faddeev equation one obtains the nucleon's mass, and also 
its Faddeev amplitude which is a valuable intuitive tool. It is noteworthy that 
even a rudimentary covariant Faddeev equation model, based on the presence of 
diquark correlations within the nucleon, yields a matrix- valued amplitude that, 
in the nucleon's rest frame, corresponds to a relativistic wave function with a 
material lower component; i.e., a wave function with "p-wave," and, indeed, a d- 
wave" correlations, too |26| . Nonzero quark orbital angular momentum in the 
nucleon is a straightforward outcome of a Poincare covariant description. 

While some issues remain unresolved j57J|55], contemporary data [2*^1 I5U1 
l3T] suggest that a single dipole mass cannot simultaneously characterise the 
Q 2 -dependence of both the proton's electric and magnetic form factors. This 
possibility was evident in the Faddeev-amplitude-based calculations of Ref . ^0J , 
as emphasised in Ref. Moreover, it can be inferred from Ref. |32| that 

this experimental result is an essential consequence of a nonperturbative and 
Poincare covariant representation of the proton as a bound state. This may be 
exemplified through the role played by pseudovector components of the pion's 
Bethe-Salpeter amplitude, which are connected with the presence of quark orbital 
angular momentum in the pion. These pseudovector amplitudes are necessarily 
nonzero [3S] and responsible for the large-Q 2 behaviour of the electromagnetic 
pion form factor |341 I35| . 

Herein we calculate the proton's elastic electromagnetic form factors, using a 
product Ansatz for the proton's Faddeev amplitude and a spectator approxima- 
tion to describe elastic electromagnetic scattering from the nucleon. We describe 
the model in Sect. UJ and present and discuss the results in Sect. |5] Section 2] 
is an epilogue. The study furnishes a means by which we may explore and illus- 
trate the points outlined above. It will become apparent that existing precision 
data on the form factor ratios F| (q 2 ) / Ff (q 2 ) and G p E (q 2 )/G p M (q 2 ) are a sensitive 
probe of nonperturbative aspects of the proton's structure. 



2 Model Elements 

2.1 Dressed quarks 

There are three primary elements of our model and to begin with its specification 
we note that quarks within bound states are described by a dressed propagator 1 

S(p) = -i-f-pa v (p 2 ) +(7 S {p 2 ) = - * —, 2 = Z<J) } /T( 2 \ ■ C 1 ) 

Z7 • p A{p^) + -d(j? J i^-p + M(p z ) 



We employ a Euclidean metric wherewith the scalar product of two four vectors is 
a ■ b = X^t=i a ibi, and Hermitian Dirac-7 matrices that obey {7 M ,7i,} = 2 5 M „. 
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It is a longstanding, model-independent DSE prediction that the wave function 
renormalisation Z(p 2 ) and mass function M(p 2 ) exhibit significant momentum 
dependence for p 2 < 2GeV 2 whose origin is nonperturbative [UElini- This be- 
haviour was recently verified in numerical simulations of quenched QCD |36| . 
and the connection between this and the full theory is analysed in Ref. |37j . 

The mass function is enhanced at infrared momenta, a feature that is an 
essential consequence of the dynamical chiral symmetry breaking (DCSB) mech- 
anism. It is also the origin of the constituent-quark mass. With increasing space- 
like p 2 , on the other hand, the mass function evolves to reproduce the asymptotic 
behaviour familiar from perturbative analyses and that behaviour is manifest for 
p 2 > lOGeV 2 [HE!- 

While numerical solutions for the dressed-quark propagator are readily ob- 
tained from a model of QCD's gap equation, the utility of an algebraic form for 
S(p) when calculations require the evaluation of numerous multidimensional in- 
tegrals is self-evident. An efficacious parametrisation, which exhibits the features 
described above and has been used extensively [UIE1 Ej- is expressed via 

as(x) = 2m F{2(x + m 2 )) + Ffax) F{b z x) [b + b 2 F(ex)] , (2) 

a v (x) = [1 - T(2{x + m 2 ))] , (3) 

x + m z L J 

with x = p 2 /X 2 , m = m/X, 

1 - e~ x 

m = , ( 4 ) 

as(x) = Xas{p 2 ) and ay(x) = X 2 (jy(p 2 ). The mass-scale, A = 0.566 GeV, and 
parameter values 2 



m bo b\ 



0.00897 0.131 2.90 0.603 0.185 



(5) 



were fixed in a least-squares fit to light-meson observables |39j . The dimensionless 
u = d current-quark mass in Eq. (JSJ) corresponds to 

m = 5.1MeV. (6) 

The parametrisation yields a Euclidean constituent-quark mass 

M u,d = 0.33 GeV, (7) 

defined as the solution of p 2 = M 2 (p 2 ) [3S1, whose magnitude is typical of that 
employed in constituent-quark models [1U]. This is an expression of DCSB, as is 
the vacuum quark condensate 

- <^>o GeV2 = ^ 3 T~2 FT - ^ n ^ = (0-221 GeV) 3 , (8) 

47r fe l 03 ^qcd 

Aqcb = 0.2 GeV. The condensate is calculated directly from its gauge invariant 
definition i33^ after making allowance for the fact that Eqs. @, © yield a chiral- 
limit quark mass function with anomalous dimension 7 m = 1. This omission of 
the additional ln(p 2 /ylQ CD )-suppression that is characteristic of QCD is merely 
a practical simplification. 

2 e = 10 -4 in Eq. (|5J acts only to decouple the large- and intermediate-p 2 domains. 
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2.2 Product Ansatz for the Faddeev amplitude 

We represent the proton as a composite of a dressed-quark and nonpointlike, 
Lorentz-scalar quark-quark correlation (diquark), and exhibit this via a product 
Ansatz for the Faddeev amplitude 

fff (»,<*, tj) = [r a \\ P[l2 y,K)]^ a2 A a \K)[S(t,P)u(P)]Z, (9) 

wherein (pi,cti,Ti) are the momentum, spin and isospin labels of the quarks 
constituting the nucleon; the spinor satisfies 

u{P) (17 • P + M) = = («7 • P + M) u(P) , (10) 

with P = pi + y»2 + P3 the nucleon's total momentum, and it is also a spinor 
in isospin space with (p + = col(l,0) for the proton and cp— = col(0, 1) for the 
neutron; and K = p 1 + p 2 =: p {12 }, P[u] = Pi ~P2, t ■= {~P{12} + 2p3)/3. 

In Eq. Q, A 0+ (K) is a pseudoparticle propagator for the scalar diquark 
formed from quarks 1 and 2, and r 0+ is a Bethe-Salpeter-like amplitude describ- 
ing their relative momentum correlation. These functions can be obtained from 
an analysis of the quark-quark scattering matrix, as explained in Ref. [25] . How- 
ever, we have already chosen to simplify our calculations by parametrising S(p), 
and hence we follow Refs. jlOl 111] and also employ that expedient herein, using 

A 0+ (K) = -l-F^/ufa), (11) 

r 0+ (k;K) = j^HCijs^Fik 2 /^), (12) 

with T defined in Eq. (JIJ, C = 7274 the charge conjugation matrix, t 2 the 2x2 
Pauli isospin matrix, and (H C3 ) C1C2 = e clC2C3 , 01,2,3 = 1)2,3, describing the com- 
pletely antisymmetric colour structure of a 3 diquark. 3 The parameters are a 
width, u! +, and a pseudoparticle mass, m +, which have ready physical inter- 
pretations: the length l + = 1/uj + is a measure of the mean separation between 
the quarks in the scalar diquark; and the distance Aq+ = l/m + represents the 
range over which a diquark correlation in this channel can persist inside the 
nucleon. 

With the elements described hitherto it is possible to derive a Poincare co- 
variant Faddeev equation whose solution yields S, a 4 x 4 Dirac matrix that 
describes the relative quark-diquark momentum correlation. The complete nu- 
cleon amplitude then follows: 

& = ^ + ^ + > (13) 

where the subscript identifies the dormant quark and, e.g., are obtained from 
•F3 by a uniform, cyclic permutation of all the quark labels. The general form of 

3 In Eq. 112|l . A/" 0+ is a calculated normalisation constant which ensures that a (wd)-diquark has 
electric charge fraction (1/3) for K 2 — -irij + . TV*, to appear in Eq. 1)15^. is analogous: it is 
the calculated normalisation constant that ensures the proton has unit charge. 
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S is discussed at length in Ref. [2E1 with the conclusion that the positive energy 
solution can be written 

S{£; P) = f 1 (l;P)I D + ± ! (i 1 .£-t.p J D ) f 2 (£; P) , (14) 

where (Ij^) rs = S rs , P 2 = —1. In the nucleon's rest frame, /12, respectively, 
describe the upper, lower component of the spinor amplitude S{£; P) u(P). 

Again, calculations are simplified if one employs an algebraic parametrisation 
of S, and the form jlj 



1 



S(£; P) = Tr T{£ 2 ju? M ) Id-^{h-£-£-PId) (15) 



R 



M 



is efficacious. In writing this one exploits results of the Faddeev equation calcula- 
tions reported in Refs. |25| I26j . which establish the fidelity of the approximations 
h(£;P) « Rfi(i\P) and fi(£;P) « fi(£ 2 ]P 2 ). The Ansatz involves two param- 
eters: a width oJ q {qq} and ratio R. The former can be associated with a length- 
scale Iqiqqi = l/u q { qq y, which measures the quark-diquark separation. The latter 
gauges the importance of the lower component of the positive energy nucleon's 
spinor amplitude. Its magnitude increases with increasing R. (The strength of 
the lower component of the nucleon's Faddeev wave function is determined by 
R but does not vanish for R = 0.) In realistic Faddeev equation studies of the 
nucleon l q i qq \ > Z + /2 ~ 0.2 fm and R ~ 0.5 



2.3 Dressed-quark-photon coupling 

A calculation of the electromagnetic interaction of a composite particle cannot 
proceed without an understanding of the coupling between the photon and the 
bound state's constituents. If those constituents are dressed then the coupling is 
not pointlike. Indeed, it is readily apparent that with quarks dressed as described 
in Sec. 12.11 only a dressed-quark-photon vertex, r^, can satisfy the vector Ward- 
Takahashi identity: 

q^P^xM = S~\£i) - S-\£ 2 ) , (16) 

where q = £\ — £2 is the photon momentum flowing into the vertex. This is 
illustrated with particular emphasis in Refs. [S3 113 ESI EH S3 EE1 E3| , which 
consider effects associated with the Abelian anomaly. The Ward-Takahashi iden- 
tity is only one of many constraints that apply to i - ^ in a renormalisable quantum 
field theory and these have been explored extensively in Refs. |48| I49j . 

The dressed-quark-photon vertex, a three-point Schwinger function, can be 
obtained by solving an inhomogeneous Bethe-Salpeter equation. This was the 
procedure adopted in the DSE calculation [7j that successfully predicted the 
electromagnetic pion form factor (HOI- For our purposes, however, it is enough to 
follow Ref. 39 and employ the algebraic parametrisation jSJ 



ir iM {£ 1 ,£ 2 ) = iZA(£\A) 7m + (h + h)^ [i»7 ' (h + h) A A (£ 2 ,£ 2 ) + A B (£ 2 , £ 2 2 )\ ; 

(17) 
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with 



£ F (li£l) = l[F(£l)+F(ll)), A F {l\A) 



i\-t\ 



(18) 



where F = A,B; i.e., the scalar functions in Eq. It is critical that r"„ in 
Eq. 1)17(1 satisfies Eq. (|16|) and very useful that it is completely determined by 
the dressed-quark propagator. Improvements to this Ansatz modify results by 
< 10%, as illustrated, e.g., in Refs. jHHESl- 

Equation (|17|) entails that dressed-quarks do not respond as point particles 
to low momentum transfer probes. This observation qualitatively supports an 
assumption employed in some relativistic constituent quark models ^3 E3J • 
An unambiguous quantitative connection is difficult because the definition of 
constituent-quark degrees of freedom depends on a model's formulation. It may 
nevertheless be worth noting that quark dressing disappears with increasing 
spacelike q 2 in QCD. Hence, in the ultraviolet, the dressed-quark's Dirac form 
factor must approach one (with only In q 2 corrections) and its Pauli form factor 
must vanish; viz., the interaction becomes pointlike in this limit. With the pa- 
rameter values in Eq. ©, this evolution of the Dirac and Pauli form factors may 
be characterised by monopole ranges r\ ~ 0.25 fm and r<i ~ 0.35 fm, respectively. 

2.4 Commentary 

We have completely specified a covariant model of the nucleon as a bound state of 
a dressed-quark and nonpointlike scalar quark-quark correlation. This algebraic 
Ansatz has four parameters: m + and co + introduced in Eqs. (|TT|) . ((T2*|) to charac- 
terise the diquark; and ^ q { qq } and R in Eq. (|T5|) . which express prominent features 
of the nucleon's spinor. The dressed-quark propagator and dressed-quark-photon 
vertex are fixed. 

In contemplating such a model one may ask whether it can supply an accurate 
description of the nucleon's electromagnetic form factors. A priori, the answer 
is unknown but it is supplied by straightforward calculations. 

A more important question, however, is whether the model should be accu- 
rate. In this case the answer is no. Reference j2S! emphasises that no picture of 
the nucleon is veracious if it neglects axial-vector diquark correlations and the 
nucleon's pion "cloud." Thus our simple model must be incomplete. Fortunately, 
estimates exist of the contributions made by these terms to the nucleon's electro- 
magnetic properties |12 | I55 [ 156 } I57j . and in the following they are used to inform 
the model's application. 

3 Calculated Form Factors 

The nucleon's electromagnetic current is 



where q = P' — P and A^ is the nucleon-photon vertex described in the appendix. 
In Eq. (J20[). F\ and F<i are, respectively, the Dirac and Pauli electromagnetic form 




ieu(P') 1 ^F l {q 2 ) + — cv q v F 2 (q 2 ) u(P) 



(19) 



(20) 



J. C. R. Bloch, et al. 



7 



Table 1. Fitted model parameters and calculated nucleon static properties. Experi- 
mental values are provided for comparison in the last row. NB. The values indicate: 
A 0+ = l/m 0+ ~ ifm, ; + = l/w 0+ ~ ifm, l q{qq} = l/Lu q{qq} ~ | fm. 





Parameters 




Calculated Static Properties 




R 


m 0+ (GeV) uj 0+ 


^q{qq} 


(rl) 2 (fm 2 ) 






Hn 


0.25 


0.75 1.50 


0.33 


(0.65) 2 


-(0.38) 2 


2.58 


-1.39 


0.50 


0.77 1.42 


0.29 


(0.61) 2 


-(0.37) 2 


2.52 


-1.37 


Obs. 






(0.87) 2 


-(0.34) 2 


2.79 


-1.91 



factors. They are the primary calculated quantities, from which one obtains the 
nucleon's electric and magnetic form factors 

G E (q 2 ) = Ftf) - ^F 2 (q 2 ) , G M (q 2 ) = F^q 2 ) + F 2 (q 2 ) . (21) 

To proceed with the illustration we select two values of R, namely, R = 0.25, 
0.50. This choice is motivated by Faddeev equation studies, in which the smaller 
value is obtained by calculations that retain scalar and axial-vector diquark cor- 
relations but neglect the pion cloud, while the latter is obtained if the estimated 
effect of that cloud is incorporated 25j. Then, in each case, we fix the remaining 
three parameters by requiring a least-squares fit to 

G P E (q 2 ) = l/(l+q 2 /m 2 ) 2 , (22) 

with rrid = 1.1 GeV. This value of the dipole mass corresponds to a proton 
charge radius r p = 0.62 fm; i.e., ~ 30% smaller than the experimental value, and 
therefore leaves room deliberately for additional contributions to the nucleon's 
electromagnetic structure from axial-vector diquark correlations and a meson 
"cloud." While the leading nonanalytic contribution to r p can alone repair the 
discrepancy [33 EZ] , that does not alter the quiddity of our scheme because the 
scale of the effect is clear and redistributing strength between complementary 
contributions can be achieved merely by fine-tuning the model's parameters. 

In Table H we list the parameter values produced by this fitting procedure 
along with the calculated values of proton and neutron static properties. The 
proton's charge radius is precisely that for which we aimed. However, the re- 
maining values point to the deficiencies of a model that retains only a scalar 
diquark correlation. They confirm that in this case one is unable to obtain a 
quantitatively accurate picture of the nucleon. This is reassuring because, in 
contributing to nucleon observables, axial-vector diquark correlations primarily 
interfere constructively with the pion cloud; e.g., they both provide additional 
binding and hence act to reduce the nucleon's mass and they both act to 
increase |//j,, n |> \ r p, n \ 13 ESI E3| • Hence a model that ignores these contri- 
butions but succeeds in fitting experimental data is likely to possess spurious 
degrees of freedom. 4 The improvements necessary to make the model more real- 

4 The implementation of current conservation in the one-body current described in the appendix 




Figure 1. Left panel - Calculated proton Dirac form factor; Right panel - Calculated proton 
Pauli form factor. Solid line, R = 0.25; dashed, R = 0.50. The associated model parameters are 
listed in Table □ 



istic are therefore plain. In their absence it is nevertheless possible to illustrate 
important points. 

A first observation relates to what may be called the nucleon's "quark core." 
In an internally consistent model it is always possible to identify the relative 
strength of various contributions to a physical observable. For example, Ref. 
|25| employs a rainbow-dressed quark and ladder-bound meson basis, within 
which the quark core contributes approximately 85% of the nucleon's mass. In 
the present case an estimate of the core's spatial extent is afforded by l q { qq } = 
~ | fm, which is commensurate with our calculated core contribution 
to the proton's charge radius, Table ^ It will also be evident that this scale is 
consistent with estimates of the magnitude of meson-loop contributions to the 
proton's charge radius determined from lattice-QCD simulations |551 157j . 

We depict the proton's Dirac form factor in the left panel of Fig. ^ wherein 
it is clear that the shift in parameter values has little observable impact. That 



is too simple to allow a fair description of G\. In this case it misses important cancellations 
and hence we report anomalously large values for \r„\. A realistic description requires a more 
complex current which includes fully- fledged seagull terms |12) . The simple current is adequate 
for the remaining form factors because such destructive interference is either absent or markedly 
less important |10|. 
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q 2 (GeV) 2 

Figure 2. Calculated proton electric form factor: solid line, R = 0.25; dot-dashed, R = 0.50 
- the associated model parameters are listed in Table □ For comparison, the lighter dashed 
line is G p E (q 2 ) = 1/(1 + 5 2 / m Lp-p) 2 i wi cmp _ p = 0.84 GeV; viz., a dipole fit to available proton 
data. 



is also true for the Pauli form factor, which is plotted in the right panel. These 
two functions together give the proton's electric form factor, via Eq. (|21|) . and 
that is shown in Fig. "2J 

We plot the calculated ratio: Fz/kFi, k = fj, p — 1, in Fig. [""J along with 
modern experimental data |29| \M\. In addition, we draw the calculated result 
for this ratio obtained using the parameter values of the scalar diquark Ansatz 
employed in Ref. (THj. These values: R = 0.0, and (in GeV) m + = 0.63, uj + = 
1.4, uj q { qq y = 0.2, were fixed via a least-squares fit to Eq. ((2*2*)) but with = 
0.84 GeV, which gives r p = 0.87 fm. The figure illustrates that the ratio decreases 
smoothly with increasing R, bracketing the data, and thereby suggests that even 
the rudimentary Ansatz is capable of accurately describing the data. Indeed, 
the R = 0.25 parameter set might even be considered a good representation. 
However, this is not the conclusion we draw. Rather, the result demonstrates 
that the available data on this ratio is sensitive to model-dependent details. 

The ratio ix p G p E /G p M is depicted in Fig. [Hand therein the last statement is 
amplified. The proton's electric form factor is a difference, Eq. (|21[) . and that 
accentuates its sensitivity. On the domain for which data is available, this ratio is 
particularly responsive to model details, a result also conspicuous in Ref. [12]. It 
is plain that three parameter sets, which are reasonable and differ modestly when 
compared through the ratio F^/kFi, appear vastly different in the comparison 
presented in Fig. [3] Moreover, because of continuity, it is clear that one could 
tune the model parameters to fit the data on this ratio. However, what might 
be considered success in that endeavour could easily be achieved through results 
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1 2 3 4 5 6 

q 2 (GeV 2 ) 



Figure 3. Calculated ratio F$/(KFf), k = (jl p — 1: solid line, R = 0.25; dashed, R = 0.50 - 
the associated model parameters are listed in Table Q The short-dashed line was obtained with 
the model parameters in Ref. |10|: R = 0.0, and (in GeV) m + = 0.63, lo + = 1.0, u q { qq } = 0.2. 
Data: boxes, Ref.|29_j; circles, Ref.|31j. 

for G P E and G P M individually which both disagree with the data. 
4 Summary and Discussion 

We calculated the proton's elastic electromagnetic form factors using a rudimen- 
tary Ansatz for the nucleon's Poincare covariant Faddeev amplitude that repre- 
sents the proton as a composite of a confined quark and confined nonpointlike 
scalar diquark. All such models give a Faddeev wave function that corresponds 
to a nucleon spinor with a sizeable lower component in the rest frame. 

This study indicates that on the domain of q 2 accessible in modern precision 
experiments these form factors are a sensitive probe of nonperturbative dynam- 
ics. The calculated pointwise forms express a dependence on the length-scales 
that characterise nonperturbative phenomena, such as bound state extent, dress- 
ing of quark and gluon propagators, and meson cloud effects. This is precisely 
analogous to the current status of the pion's electromagnetic form factor, for 
which the behaviour predicted in a straightforward application of perturbative 
QCD is not unambiguously evident until q 2 > 15 GeV 2 

The ratio fi p G p E /G p M is particularly sensitive to infrared dynamics because, 
as q 2 increases, G E (q 2 ) is a difference of small quantities. However, this ratio 
should not be considered in isolation because it is possible to reproduce the ex- 
perimental data using a model that simultaneously provides a poor description 
of the individual form factors. That is also true of the ratio F p / kF p but this 
combination is less responsive to model particulars because the Dirac and Pauli 



J. C. R. Bloch, et al. 



11 




Figure 4. Calculated ratio /i p G p E (q 2 ) /G p M (q 2 ): solid line, R = 0.25; dashed, R = 0.50 - the 
associated model parameters are listed in Table Q The lighter short-dashed line was obtained 
with the model parameters in Ref. |10|: R = 0.0, and (in GeV) m + = 0.63, u + = 1.0, 
LO q { qq } = 0.2. Data: boxes, Ref . |29j : diamonds, Ref.|30): circles, R,ef.|31). 



form factors are positive and fall uniformly to zero, with a momentum depen- 
dence at asymptotically large momenta given by perturbative QCD [S^ll^U]. In 
the absence of a veracious theoretical understanding of the nucleon, we view the 
latter ratio as a more sensible constraint on contemporary studies. 

It is apparent that much can be learnt about long-range dynamics in QCD 
from existing and forthcoming accurate data on nucleon form factors. In con- 
straining systematic QCD-based calculations, one can hope, for example, to see 
an evolution from the domain on which meson cloud effects are important to that 
whereupon observables are dominated by the confined quark core. This could be 
elucidated by improving the present study; viz., basing it on a solution of the 
Poincare covariant Faddeev equation with axial-vector diquark correlations, in- 
stead of using an Ansatz for the Faddeev amplitude, and explicitly including 
meson cloud contributions. The form factors would then be tied directly to as- 
sumptions about the nature of quark-gluon dynamics in the nucleon. 
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Figure A.l. One-body current obtained from the product Ansatz of Eq. JJ. Solid line, dressed- 
quark propagator S(p); dashed line, diquark propagator A {>+ (K). The amputated nucleon ver- 
tex is ^3 + in every case. 



Appendix A: Nucleon-Photon Vertex 

We use an Ansatz for 3*3 in the nucleon's Faddeev amplitude, Eq. 1131 . from which a properly 
antisymmetrised one-body vertex can be constructed via the method outlined in Ref. 



(A.l) 



wherein: 

Al(q, P) = 3 S(£ - |«; P) A 0+ (K) S(£; P) A^P3 + q, P3 ) , (A.2) 

with K = I + ~P, P3 = ip - t, A*(<i,«i) = 8{£i) r?{h,£2) S(£ 2 ); and 

Al(q,P) = 6 y^^I(^yi ^(Pi +Q,P2,P3) 0( Pl ,p 2 ,p 3 )tT DF [^(pi +g,Pi)S(p 2 )] 5(p 3 ) , 

(A.3) 

1 1 

where pi = — it + fc, P2 = — if — k, 6 = £ C ic2c 3 £cic2C3 is the colour contraction, and 

M(pi,p2,ps) = ^ 0+ (P{i2})P +(^P[i2];P{i2})5(i[p {12 } -2p 3 ];P). (A.4) 

A? describes the photon probing the structure of the scalar diquark correlation, and contributes 
equally to both the proton and neutron. The remaining terms are 

i 4 fc dH 

(2tt) 4 (2tt) 4 

x 0( Pl + q,p 3 ,p 2 ) S(p 2 ) {iT 2 y A®( Pl , Pl + q) {ir 2 ) 5(p 3 ) , (A.5) 
where "t" denotes matrix transpose, and 

/d^ k d^ £ 
7^4 ^(Pi.fa-Pa + ?) ^(pi,P3,P3) ^(pa + 9,P2) S(pi) S(p 3 ) , 
(27Tj^ (27T)* 



l3 , „. „ f d A k d 4 £ n , 
^(g,P) = 6 y ^(pr,p 2 ,P3) 
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(A.6) 

J2%Y (2tt)* Q(j>1,P3 + Q'PtW&i-'PZ'P 3 ) S (P 2 ~) S ^ A ^P 3 + 9'^ 3 ) ' 

(A.7) 

We illustrate these five terms in Fig. IA.1I They are in one-to-one correspondence with those 
considered in Ref. |58) . with the bottom two diagrams, representing A^' 5 , being progenitors of 
the "seagull" terms exploited therein to ensure current conservation. 

Our results are obtained by evaluating these integrals using Monte-Carlo methods and the 
input specified in Eqs. ©, ®, GD, E3 and CU- 
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